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We study the effect of the lattice structure on the spin-fluctuation mediated superconductivity 
in the iron pnictides adopting the five-band models of several virtual lattice structures of LaFeAsO 
as well as actual materials such as NdFeAsO and LaFePO obtained from the maximally-localized 
Wannier orbitals. Random phase approximation is applied to the models to solve the Eliashberg 
equation. This reveals that the gap function and the strength of the superconducting instability 
are determined by the cooperation or competition among multiple spin fluctuation modes arising 
from several nestings among disconnected pieces of the Fermi surface, which is affected by the 
lattice structure. Specifically, the appearance of the Fermi surface 7 around (vr, tt) in the unfolded 
Brillouin zone is sensitive to the pnictogen height ftpn measured from the Fe plane, where /ipn is 
shown to act as a switch between high- Tc nodeless and low- Tc nodal pairings. We also find that 
reduction in the lattice constants generally suppresses superconductivity. We can then combine 
these to obtain a generic superconducting phase diagram against the pnictogen height and lattice 
constant. This suggests that NdFeAsO is expected to exhibit a fully-gapped, sign-reversing s-wave 
superconductivity with a higher Tc than in LaFeAsO, while a nodal pairing with a low Tc is expected 
for LaFePO, which is consistent with experiments. 
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I. INTRODUCTION 

The discovery of superconductivity in the iron-based 
compounds by Hosono's group^ and subsequent increase 
in the transition temperature (Tc) exceeding 50 K— in the 
same family of compounds are seminal not only because 
of high values of Tc, but also because this poses a fun- 
damental question on electronic mechanisms of high Tc 
superconductivity in a wider class of compounds other 
than cuprates. 

Theoretically, a phonon mechanism was shown to be 
unlikely for this system,"^ and a spin-fluctuation mediated 
pairing has been proposed from the very early stage of 
the studyi^i^i^ In these studies the nesting between dis- 
connected pieces (pockets) of the Fermi surface is shown 
to induce spin fluctuations associated with the nesting 
vector. This can give rise to a superconducting gap, 
which is basically s-wave but changes sign between dif- 
ferent pockets, hence termed as sibwave or sign-reversing 
s-wave first proposed by Mazin et alA{see Fig|4]). Al- 
though recent experimental as well as theoretical stud- 
ies suggest that the magnetism in the undoped material 
is not driven solely by Fermi surface nesting,^ the spin 
fluctuation originating from the nesting has been consid- 
ered to be a possible origin of the pairing interaction by 
a number of authors^i^^^^ii^ii^^i^^ii^ 



In particular, the present authors with Tanaka and 
Kontani have constructed a minimal model, which has 
turned out to be five-band, for LaFeAsO based on first 
principles calculation, and investigated spin-fluctuation 
mediated superconductivity with random-phase approx- 
imation (RPA),^ii^i^i2i In that study it was pointed out 
that, along with the sign-reversing s-wave, a c?-wave pair- 
ing can also be a candidate depending on the band flUing. 
Our five-band model has subsequently been adopted in 
various studies, among which are a perturbation study 
by Nomura^^, a fluctuation exchange study by Ikeda— , 
and a functional renormalization group study by Wang 
et ali^ An analysis on the normal state spin dynam- 
ics has also been performed using our five-band modelj^ 
where good agreement with inelastic neutron scattering 
experiments^!^ has been obtained. On the other hand, 
Graser et al^'^ recently applied RPA to a five-band 
model that is based on a band structure obtained by Cao 
et al2-^ The study finds that a sign-reversing s-wave that 
has nodes intersecting the Fermi surface closely competes 
with d-wave pairing (see FigU]). It has further been pro- 
posed in ref , 28 that this nodal s-wave pairing is intrinsic 
to the iron pnictide superconductors, while a full gap can 
occur as a consequence of the presence of impurities. It 
is also worth noting that a competition or mixture of 
sign-reversing s and d pairings have also been discussed 
on the basis of a 16 band d-p model^^ and a two-orbital 
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exchange coupling model (J1-J2 model) rSE It is the pur- 
pose of the present paper to explore systematically the 
material- and structure-dependences on the strength and 
gap symmetry of superconductivity in terms of the five- 
band model. 

Experimentally, the fully-gapped, sign-reversing s- 
wave scenario is consistent with a number of mea- 
surements on arsenides, such as angle-resolved pho- 
toemission spectroscopy (ARPES)r^i^ penetration 
depth measurements,'^^ and muon spin relaxation 
(/^SR ) "S i ^.4 i i?5 | i?^., which suggest that the gap is open on 
the whole Fermi surface, although the magnitude of the 
gap may vary along the surface. The fully-gapped, sign- 
reversing s is also consistent with some neutron scatter- 
ing results^i^^ that observe a resonance peak predicted 
theoreticallyi^2ii2iiii On the other hand, the weak effect 
of nonmagnetic impurities such as Co on or even the 
appearance of superconductivity upon Co doping^^^i^ 
has cast doubt on the sign-reversing gap, but some theo- 
retical studies^i^i^ have shown that these experiments 
in fact do not necessarily contradict with the s±. In 
particular, Senga and Kontani showed that the effect of 
the inter-pocket scattering due to nonmagnetic impuri- 
ties becomes irrelevant in the unitarity limit. 

One interesting and important feature in the iron pnic- 
tides is the unusually strong dependence of Tc on ma- 
terials, which ranges from ~ 5K in LaFePO''^ to 55K 
in SmPeAsO^ even within the same group of elements. 
More systematically, Lee et al^ pointed out that we 
can parametrize the value of Tc in terms of the Fe-Pn-Fe 
(Pn: pnictogen) bond angle a, where Tc seems to have a 
peak around the bond angle (~ 109) at which the pnic- 
togens form a regular tetrahedron, while Tc is low for 
materials with large a such as LaFePO. The importance 
of the bond angle has also been pointed out by Zhao 
et al.^^ On the other hand, it has also been shown by 
Miyazawa et al. that the chemical trend for LnFeAsO 
has the maximum Tc increasing with the decreasing lat- 
tice constant a for Ln=La — > Nd, but nearly constant 
for Ln=Nd — > Dyi^ Pressure effects have also been ex- 
perimentally elaborated. For LaFeAsO, Tc first increases 
with pressure, but then decreases when the pressure be- 
comes too largCf^i^i^i^ which is contrasted with mate- 
rials having Tc > 50 K at ambient pressure such as Nd- 
FeAsO for which Tc monotonically and rapidly decreases 
with pressure. All these experimental results indicate 
that Tc is unusually sensitive to the lattice structure in 
the iron-based compounds. 

If we move on to the symmetry of the superconduct- 
ing gap, we have various pieces of experimental evidence 
for strong material dependence as well: while a num- 
ber of experiments on arsenides suggests that the gap 
is fully open on the Fermi surface as mentioned, a re- 
cent penetration depth measurement on LaFePO shows 
that there are nodes in the superconducting gap i^^i^^ 
Arsenides and LaFePO also exhibit sharp contrast in 
nuclear magnetic resonance (NMR) experiments. In 
LaFeAsO, some experiments show that the NMR re- 



laxation rate 1/Ti has no coherence peaks, and decays 
as below T^^i^i^iSi^i^ while a recent experiment 
by Kobayashi et al^ indicates a more rapid decay. In 
LaFePO, 1/Ti below Tc is shown to decay even more 
slowly than above Tc"^ These results strongly suggest 
that even the symmetry of the superconducting gap can 
be unusually sensitive to the lattice structure. 

So the crucial theoretical question is how we can under- 
stand these sensitivities. In analyzing the structure de- 
pendence, there is one factor to which we have to pay at- 
tention. Previous theoretical studies have shown that the 
position of the pnictogen with respect to the Fe plane af- 
fects the band structure, in particular the character of the 
bands that lie close to the Fermi level near the F point (in 
the folded Brillouin zone) as well as the band width (see 
FigO^i^il^ Local spin-density approximation studies 
have shown that the tendency towards magnetism be- 
comes stronger when the pnictogen lies farther from the 
Fe plane, which is expected to enhance superconductivity 
if the pairing is mediated by spin fluctuationsj^ 

Given this background, in the present study we in- 
vestigate the lattice structure dependence of the spin- 
fluctuation mediated superconductivity, where we con- 
struct five-band models for several virtual lattice struc- 
tures of LaFeAsO as well as actual materials such as Nd- 
FeAsO and LaFePO, and apply RPA to solve the Eliash- 
berg equation. We shall show that the position of the 
pnictogen is indeed the key factor that determines both 
Tc and the form of the superconducting gap, namely, the 
"pnictogen height" above the Fe plane (Fig[T]) can act as 
a switch between a high-Tc, fully- gapped, sign-reversing 
s-wave and a low-Tc, nodal (s- or d-wave) pairings. We 
also show that the reduction in the lattice constant is 
generally unfavorable for superconductivity. Combining 
these results for the lattice structure dependence, we then 
obtain a generic "phase diagram" against the pnictogen 
height and the lattice constants. Based on the phase di- 
agram, we argue that the systematic dependence of Tc 
against the bond angle found in refj5(l can be accounted 
for as a combined effect of the pnictogen height and the 
lattice constants. In order to get higher Tc, we propose 
to seek for materials that have high position of the pnic- 
togen and large lattice constants simultaneously. 




FIG. 1: (Color online) Lattice structure of one Fe-Pn layer, 
with the pnictogen height indicated. 
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FIG. 2: (Color online) (a) The original (dashed lines) and 
reduced (solid) unit cells with • (Fe), V (As below the Fe 
plane) and A (above Fe). (b) The band structure (left) of the 
five-band model for LaFeAsO, and the Fermi surface (right) 
aX kz = Q for n = 6.1. The main orbital characters of some 
portions of the bands and the Fermi surface are indicated. 
The dashed horizontal line in the band structure indicates the 
Fermi level for n = 6.1. The short arrow in the band structure 
indicates the position of the Dirac cone closest to the Fermi 
level. The gray areas in the Fermi surface around the zone 
corners represent the 7 Fermi surface, (c) The portion of the 
band that has mainly the dx2_Y2 orbital character. 



then construct the maxinially-localize(i Wannier func- 
tions (MLWFs)^3^ These MLWFs, centered at the two 
Fe sites in the unit cell, have five orbital symmetries 
(dsz^-R-^, dxz, dyz, (ix2-y2, und dxY, where X,Y,Z 
refer to those for the unit cell with two Fe sites as shown 
in Figl^Ja)). The two Wannier orbitals in each unit cell 
are equivalent in that each Fe atom has the same local 
arrangement of surrounding atoms. We can then take 
a unit cell that contains only one orbital (for each or- 
bital symmetry) by unfolding the Brillouin zone, and we 
end up with an effective five-band model on a square lat- 
tice, where x and y axes are rotated by 45 degrees from 
X-Y , to which we refer for all the wave vectors here- 
after. We define the band filling n as the number of 
electrons/number of sites (e.g., n = 10 for a full filling). 
The doping level x in LaFeAsO i-ajF^, is related to the 
band filling as n — Q -\- x. 

The five bands are heavily entangled as shown in 
FiglUb) refiecting the strong hybridization of the five 2>d 
orbitals, which is physically due to the tetrahedral coor- 
dination of As atoms around Fe. Hence we conclude that 
the minimal electronic model requires all the five bandsJ^ 
In Fig[2lb), the Fermi surface at fcz = for n = 6.1 (cor- 
responding to X ~ 0.1) is shown in the unfolded Bril- 
louin zone. The Fermi surface consists of four pieces 
(pockets in 2D): two concentric hole pockets (denoted 
as ai, 0:2) around (kx,ky) = (0,0), two electron pockets 
around (tt, 0) (/3i) or (0,7r) [(32), respectively. Besides 
these pieces of the Fermi surface, there is a portion of the 
band around (tt, tt) that touches the Ep at n = 6.1, so 
that this portion acts as a "quasi Fermi surface" (which 
we call 7). As for the orbital character, a and some por- 
tions of [3 near the Brillouin zone edge have mainly dxz 
and dyz character, while the portions of /3 away from the 
Brillouin zone edge and 7 have mainly dx^-v'^ orbital 
character. An interesting feature in the band structure 
is the presence of Dirac cones, i.e., places where the up- 
per and the lower bands make a conical contact i^i^i The 
ones closest to the Fermi level correspond to the crossing 
points of the dx^-y'^ the dxz /dyz bands below the 
[3 Fermi surface. 



III. MANY-BODY HAMILTONIAN AND 
RANDOM-PHASE APPROXIMATION 



II. THE BAND STRUCTURE AND THE FERMI 
SURFACE 



LaFeAsO has a tetragonal layered structure, where Fe 
atoms form a square lattice in each layer, which is sand- 
wiched by As atoms (Figs lll2r a)). Due to the tetra- 
hedral coordination of As, there are two Fe atoms per 
unit cell. The experimentally determined lattice con- 
stants are a = 4.036A and c = 8.739A, with two internal 
coordinates ZLa = 0.142 and zas = 0.6512ii We have 
obtained the band structure (Figl^I^b) ) with the local 
density approximation with a plane-wave basis^. We 



For the many body part of the Hamiltonian, we con- 
sider the standard interaction terms that comprise the 
intra-orbital Coulomb U , the inter-orbital Coulomb [/', 
the Hund's coupling J, and the pair- hopping J'. The 
many body Hamiltonian then reads 
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where i,j denote the sites and /x, ^ the (five d) orbitals, 
and t^!^ is the obtained in the previous section. The or- 
bitals (i3Z2__R2, dxz, dyz, dx^-Y^, and dxY are labeled 
as = 1,2,3,4, and 5, respectively. As for the electron- 
electron interactions, there have been theoretical studies 
that estimate the parameter values. Some studies give 
U = 2.2- 3.3 and J = 0.3 - O.e'^^-'^^ in units of eV, while 
others have U ^ J 7"^ Here we assume that U > J and 
take the values U = 1.2, U' = 0.9, and J = J' = 0.15. 
We also examine orbital-dependent interactions as intro- 
duced in sectior fVCl We have taken the values somewhat 
smaller than those obtained in ref lTsU?^ because the self- 
energy correction is not taken into account in the present 
RPA calculation, so that small interaction parameters are 
needed to avoid magnetic ordering at high temperatures. 

Having constructed the model, we move on to the five- 
band RPA calculation, where the modification of the 
band structure due to the self-energy correction is not 
taken into account. Multiorbital RPA is described in e.g. 
ref .78,79. . In the present case. Green's function Gi„i{k) 
[k = (k,iu)n)) is a 5 X 5 matrix. The irreducible suscep- 
tibility matrix 



Xi 



[li = 1,...,5) has 5'* components, and the spin and the 
charge (orbital) susceptibility matrices are obtained from 
matrix equations. 



Xs{q) 



Xc{q) 



X'iq) 



1 + Cx°(«) 



(3) 



(4) 



where 



L/, ^1 — ^2 — '3 — h 

U', h=h^l2= h 

J, h = h h = h 

J', h = I4 ^ h — h, 



U h — I2 — — 

-U' + J ll=h^l2= k 

2U' - J, h=h^l^ = U 

J' h ~ h h — h- 



(5) 



(6) 



We denote the largest eigenvalue of the spin (charge) sus- 
ceptibility matrix for iLUn = as Xs{k){xc{k)). 

The Green's function and the effective singlet pairing 
interaction, 
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v^q) = T;SUq)s - 7^cUq)c + t^{s + c), (?) 



are plugged into the linearized Eliashberg equation, 

X G';,i,(fc - q)<t>i.j^{k - q)Gi,i^{q - k). (8) 

The 5x5 matrix gap function (pim in the orbital represen- 
tation along with the associated eigenvalue A is obtained 
by solving this equation. The gap function can be trans- 
formed into the band representation with a unitary trans- 
formation. The calculation is performed at T = 0.02 eV 
taking a three dimensional fc-point mesh of 32 x 32 x 4 and 
512 Matsubara frequencies. All the results for the spin 
susceptibility and the superconducting gap will be pre- 
sented for the lowest Matsubara frequency and at k^ = 
or qz = 0. The eigenvalue of the Eliashberg equation A 
at the fixed temperature of 0.02eV will be adopted as 
a measure of the strength of the superconducting insta- 
bility, since directly obtaining Tc, especially for low Tc 
systems, requires more fc-point meshes and Matsubara 
frequencies. 



IV. ORBITAL-DEPENDENT NESTING AND 
THE PAIRING SYMMETRY COMPETITION 

Let us first look in Figl3]at the result for the (orbital- 
diagonal components of) spin susceptibility, Xs3333 and 
Xs4444: which are the two largest components. Xs3333 has 
peaks solely around (tt, 0) and (0, tt), which reflects the 
nesting between dxz, dyz portions of a and /3 pockets as 
shown in the lower panel of Fig[31 where the thickness of 
the Fermi surface represents the strength of the dx^-v^ 
or dxz /dyz characters. On the other hand, XsAaaa has 
peaks around (tt, 0), (0, tt) and (tt, 7r/2), (7r/2, tt). The for- 
mer is due to the nesting between the 7 pocket and the 
dx'^-y'^ portion of the (3 pocket, while the latter origi- 
nates from the nesting between the dx'^~y'^ portion of 
the /3i and/32-^'21^'26^ 

The superconducting gap should be determined by the 
cooperation or competition between the multiple nestings 
mentioned above. Specifically, the a-[3 and 7-/? nestings 
tend to favor the fully-gapped, sign-reversing s-wave, in 
which the gap changes sign between a and (3 but has 
a constant sign on each pocket as shown in Figdl^ On 
the other hand, /?i-/?2 nesting tends to change the sign 
of the gap between these pockets, which can result in 
either d-wave or an s-wave pairing with nodes on the (3 
Fermi surface, as shown schematically in FigUiiii^. For 
the band structure of LaFeAsO (obtained by using the 
experimentally determined lattice structure), the sign- 
reversing s-wave with no nodes intersecting the Fermi 
surface dominates for the present set of parameter values 
with n = 6.1 as shown in FigO^ The eigenvalue of the 
Eliashberg equation at T = 0.02eV is A = 0.90 for s-wave, 
against A = 0.54 for d-wave. 

As for the band filling dependence, we plot the eigen- 
value of the Eliashberg equation of s- and d— wave pair- 
ings in FiglSJa). We can see for the band structure 




FIG. 3: (Color online) Top panels ; Diagonal components, 
Xs3333 Xs4444i of the Spin susceptibility matrix in the or- 
bital representation (3 : YZ,A : — Y"^) for the five-band 
model of LaFeAsO with n = 6.1. Bottom: Nesting is shown 
for the Fermi surface for orbitals XZ, YZ (left) and X^ — Y^ 
(right). Here the thickness of the Fermi surface represents the 
strength of the respective orbital character. 



of LaFeAsO with the present set of interaction values 
that the sign-reversing s-wave pairing with a full gap 
for each pocket dominates for the band filling n < 6.2. 
For n > 6.3, the 7 pocket becomes less effective, and 
the (tt, 0) peak in Xsna disappears as seen in the right 
panel of FiglB]Jb) . The a pocket becomes less effective as 
well, and the (tt, 0) peak in Xs3333 becomes small. Thus 
in this region, rf-wave pairing begins to dominate, and 
the subdominant s-wave gap has nodes almost touching 
the /3 as seen in the right panel of FiglHI^c). For small 
doping levels when the 7 Fermi surface is effective and 
s-wave dominates, the magnitude of the s-wave gap has 
maxima at the positions along the (3 pocket facing the F 
point (FiglHl lower left), but when the doping increases 
to n = 6.3, the s-wave gap has minima at these points. 
The gap turns out to be nearly constant on the (3 Fermi 
surface (FiglBfc), left) for the band filling n ~ 6.2. In 
this case, the gap on a (not shown) and (3 have nearly 
the same magnitude. 

Note however that the present analysis on the band 
filling dependence does not take account of the doping 
dependence of the band structure itself, which should oc- 
cur mainly due to the change of the As position caused by 
doping. We will come back to this point in sectior TVIBl 
taking NdFeAsO as an example. 
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FIG. 4: (Color online) The fully-gapped s± wave (top panel), 
the nodal s± wave (middle) , and the d-wave gap (bottom) are 
schematically shown. The solid red (dashed blue) curves rep- 
resent positive (negative) sign of the gap. The arrows indicate 
the dominating nesting vectors. 7 Fermi surface is present 
when the pnictogen height is large. 



THE EFFECT OF THE LATTICE 
STRUCTURE 



A. Pnictogen height dependence 

We now investigate the effect of the "pnictogen height 
/ipn " , namely, the distance between a pnictogen atom 
and the Fe layer, (^Pn — 0.5) x c, where zpn is the in- 
ternal coordinate of the pnictogen atom and c the c-axis 
lattice constant. As shown in previous studieSf^i^i^S zpn 
controls the relative position of dx^-Y^{= dxy) and dz^ 
bands near the F point of the folded (original) Brillouin 
zone. In the unfolded Brillouin zone, these bands ap- 
pear near (tt, tt). In Fig[71 we show the band structure in 
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TABLE I: (Color online) Materials and lattice structures considered in the present study, and the nearest and second-nearest 
neighbor hopping integrals (in eV) in the corresponding tight-binding models. Shorthands are: La (LaFeAsO), Nd (optimally 
doped NdFeAsOi_a), Nd-p (NdFeAsOi_y under the pressure of 3.8GPa), Nd-ud (underdoped NdFeAsOi-^), and P (LaFePO). 
The cases with a — 3.95A , c — 8.40A , /ias ~ 1.38A or 1.14A correspond to virtual structures of LaFeAsO. Note that 
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FIG. 5; (Color online) The s-wave (left panels) and d-wave 
(right) gap functions for the 2nd to 4th bands from top to 
bottom in the band representation for the five-band model 
of LaFeAsO with n — 6.1. Solid lines represent the Fermi 
surface, and green dashed lines the nodes in the gap. 



the unfolded Brillouin zone for virtual lattice structure 
with ZAs — 0.658 and 0.6304 with the lattice constants 
and ZLa fixed at the original values for LaFeAsO. For the 
original ^as = 0.651, the As height is /ias = 1.32A. For 
ZAs = 0.658, /lAs increases to 1.38A, which is the same 
as in the optimally doped NdFeAsO, while zas = 0.6304 



(/lAs — 1.14A) corresponds to the height of P in LaFePO. 
We see that the dx^-y^ band that forms the 7 Fermi sur- 
face around (tt, tt) rises as Has is increased, while the dz^ 
band sinks below the Fermi level. 

The reason for these shifts in the band positions can 
be understood in terms of the hopping integrals. As Has 
increases, the nearest- neighbor hopping for the dx^-y^ 
orbital decreases as shown in Table|Tl If we approximate 
the dx2-y2 portion of the bands by 

ex2-y2(fc) = — 2tx2-Y2 [cos(fc^) -t- cos(fcy)] 

—4:tx2_y2 coa{kx) cos{ky), (9) 

where t'^2_y2 stands for the second-nearest neighbor 
hopping, the energy difference between (0, 0) and (tt, tt) is 
proportional to tx2^y2 (Figl2l^c)), so that the reduction 
in tx2-y2 acts to push up the dx2-y2 band at (tt, tt). Be- 
sides the variation in the dx2-Y2 hoppings, the increase 
in hAs results in an overall reduction in the hopping inte- 
grals of other orbitals because the effective hopping path 
Fe-^ As^ Fe becomes less effective. 

The effect of varying /ias on the spin susceptibility is 
shown in Fig[5] The dyz orbital component always has 
peaks around (tt, 0),(0, tt) reflecting the a-(3 nesting. On 
the other hand, the dx2-y2 orbital component of the spin 
susceptibility Xs4444 exhibits a strong variation with Has- 
when ft-As is large and the 7 pocket is present, Xs4444 
(FiglSfa), right) is strongly peaked at (7r,0), reflecting 
the 7-/3 nesting and also the strong electron correlation 
due to the overall reduction in the band width. However, 
as /lAs is reduced, the structure around (tt, 7r/2), which 
arises from the /3i-/32 nesting, dominates (Figl5]^b), right) 
. The effect of reducing Has resembles the effect of elec- 
tron doping, but in the case of electron doping, not only 
the effect of the 7 pocket, but also that of the a becomes 
weak, so that Xs3333 is suppressed. 

The effect on the spin susceptibility in turn affects su- 
perconductivity. When Has is large, Xs4444 and Xs3333 
spin fluctuations near (tt, 0) cooperate to mediate the 
fully-gapped, sign-reversing s-wave superconductivity. 
When /lAs is small, by contrast, the (7r,7r/2) spin fluc- 
tuations begin to favor c?-wave and nodal s-wave pair- 



7 



(a) 1.0 




0.2 



6.1 6.2 6.3 

n 




(c) 



n=62, band 4 «=6.3, band 4 




FIG. 6: (Color online) (a) Eigenvalues of the Eliashberg equa- 
tion for s-wave and d-wave, respectively, for the five-band 
model of LaFeAsO plotted against the band filling n. The 
light red (or gray) symbol for the s-wave at n = 6.3 indicates 
that the gap is nearly nodal, (b) Xs3333 and Xs4444 for n = 6.3. 
(c) The s-wave gap function for band 4 with n = 6.2(left) and 
n = 6.3(right). 
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FIG. 7: (Color online) Band structure in the five-band model 
for ZAs ~ 0.658 (a) and zas = 0.6304 (b). The lattice con- 
stants are kept at the original values of LaFeAsO. X^ — Y^ and 
denote the main characters of the bands around (tt, tt, fcz). 
Dashed lines represent the Fermi energy, and the Fermi sur- 
face at fez = for n = 6.1 is shown on the right panels. 



gapless (either d-wave or nodal s-wave) superconductiv- 
ity. 



ings. In FiglHl we plot against zas (lower scale) or against 
/lAs (upper scale) the eigenvalue of the Eliashberg equa- 
tion for the s-wave and d-wave pairings, respectively. For 
large /ias where the fully-gapped, sign-reversing s-wave 
(Figlini upper right) dominates, A is large because the 
strong spin fluctuations arising from a-/3 and 7-/3 nestings 
cooperate. This is contrasted with the case of small /ias, 
where d-wave or nodal s-wave begin to dominate (Fig llOl 
left). In this region A is small because the 7-/? nesting is 
no longer effective, or to be more precise, its remaining 
effect competes with the effect of the P1-P2 nesting. It is 
worth noting that if we adopt zas = 0.638, which is the 
value determined by theoretical structure optimization, ^ 
we have the closely competing d-wave and nodal s-wave 
pairings. This is consistent with a recent RPA calcula- 
tion by Graser et al.r^ who adopted a band structure de- 
termined by a theoretical structure optimization.— The 
message of the present analysis, then, is that the pnic- 
togen height can act as a "switch" between the fully- 
gapped, sign-reversing high- Tc s-wave and the low- Tc 



B. Lattice constant dependence 

We now turn to the effect of the lattice constants. We 
consider virtual lattice structures where one of the lattice 
constants, a or c, is varied, while the pnictogen height is 
fixed at the original value for LaFeAsO. When we reduce 
the lattice constant a, we find that the nearest-neighbor 
hopping tx2-Y2 decreases, probably because the Fe-As- 
Fe angle is reduced, which may cause a suppression of the 
effective hopping via the path Fe^ As^ Fe. Nonetheless, 
most of the other in-plane hopping integrals (including 
the ones not listed in the table) are enhanced as intu- 
itively expected. On the other hand, a reduction in the 
lattice constant c is found to mainly enhance the in-plane 
dx^^Y^ hopping (apart from the obvious enhancement of 
the hopping in the c direction) . This may be because the 
As wave function is pushed toward the Fe plane for a 
reduced layer-layer distance, and the hopping between 
dx2_y2 orbitals, which is elongated in the direction of 
the As atom positions, is enhanced by this deformation. 

The eigenvalue of the Eliashberg equation is plotted 
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FIG. 8: (Color online) Xs3333 (left panels) and X34444 (right) 
for the model with (a) zas = 0.658 and (b) zas ~ 0.6304. 
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FIG. 10: (Color online) The gap functions for the model with 
ZAs = 0.6304 (left panels) or zas = 0.658 (right). From top to 
bottom: s-wave in band 3, s-wave in band 4, d-wave in band 
3, d-wave in band 4. 



FIG. 9: (Color online) s-wave and d-wave eigenvalues of the 
Eliashberg equation plotted against zas (lower scale) or /ias 
(upper scale) for n = 6.1. The lattice constants are fixed at 
the original values for LaFeAsO. For the s-wave, the open 
(solid) circles indicate that the gap is nodeless (nodal). 



as functions of the lattice constants in FigHT] We find 
that the reduction in the lattice constants tends to sup- 
press superconductivity, which can be attributed to the 
increased hopping integrals and associated suppression 
of the electron correlation. In fact, the reduction in a 
(c) enhances the XZ, YZ {X'^ — Y'^) hopping integrals, 
which leads to suppressed Xs3333 (Xs4444) as seen from 
the comparison between Fig[3] and the lower panels of 
FigHH The effect of reduced lattice constants is small for 



the competition between s and d waves (i.e., two curves 
move roughly in parallel). We note here that, although 
an increased h^s and a decreased lattice constant a both 
lead to a reduction in the Fe-As-Fe bond angle a, they 
have opposite effects on the eigenvalue of the Eliashberg 
equation (compare Figs 151 and [TT|) . 

C. Effect of the orbital-dependent interactions 

In refs lrsH?^ it is pointed out that the interaction pa- 
rameters have significant orbital dependence. In refwSl. 
the intraorbital repulsions are U = 3.27, 2.77, 2.20, and 
3.31 (eV) for dsz^-fis, dxz/YZi dx^-y^, and dxY or- 
bitals, respectively. This variation comes from the fact 
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FIG. 11: (Color online) Upper panels : s-wave and d-wave 
eigenvalues ol the Eliashberg equation plotted against a (left) 
or c (right), where /ias is fixed at the original value for 
LaFeAsO. Lower panels: Xs3333 and Xs4444 for a = 3.88A 
with c fixed at the original value (left) and c = 8.0A with a 
fixed at the original value (right). 



that each Fe 2>d orbital hybridizes with As 4p quite dif- 
ferently. Namely, while the local basis of the five-band 
model, can be represented as a linear combination 

of atomic d!^ and orbitals (as aid\ + J2j(^ijP])' 
coefficients ai,(3ij have a strong orbital dependence. For 



example, the ratio Pij/at is large for i — X — Y but 
small for e.g. i = — R^. Therefore, if we adopt a 
common value for the interaction parameters for all five 
orbitals in the five-band model, electron correlations are 
relatively overestimated for — , because p orbitals 
are more weakly correlated than d orbitals. In order to 
avoid this problem, we should use orbital dependent in- 
teractions, where the interaction for i = X'^ — Y^ is small 
compared to others. 

So we study in this section the effect of the orbital 
dependence of the interactions, taking into account the 
orbital dependence of U\ J, J' as well. Since the self- 
energy correction is not taken into account in RPA, it is 
again necessary to reduce the interactions to avoid mag- 
netic ordering at high temperatures. Here we multiply 
all the interaction parameters in refl?^ by a factor of / = 
0.42, so the intraorbital interaction is taken to be 1.37, 
1.16, 0.92, and 1.39 (eV) for dsz^-B.'^, dxz/YZ, dx^-y^, 



and dxY orbitals, respectively. Since the dx^-Y^ orbital 
has the smallest intraorbital interaction (0.92eV in the 
present calculation as compared with 1.2eV in the calcu- 
lation for orbital- independent interactions), the effect of 
the dx^-Y^ orbital is expected to be reduced compared 
with the results obtained by using the orbital indepen- 
dent interactions. 

First, we take the model for LaFeAsO to study the 
band-filling dependence. As shown in FigfTSFa). we find 
that the s-wave pairing becomes nodal for n > 6.2, i.e., 
the s-wave becomes nodal for smaller electron doping 
compared to the case with orbital independent inter- 
actions (see FiglBfa)). Also, the nodal s-wave pairing 
still slightly dominates over d-wave even sX n ~ 6.3, at 
which, for the case of orbital-independent interactions, 
the s gives way to d. If we turn to the pnictogen-height 
dependence of the eigenvalue of the Eliashberg equation 
in Fig |12r b). the s-wave is again enhanced with the in- 
creased height, but the enhancement is smaller than in 
the case of FiglHl which can be attributed to the reduc- 
tion in the 7 (namely, the dx^-Y^) Fermi surface ef- 
fect). As for the lattice constant dependence depicted in 
FigfT^c)(d). we find that the reduction in a suppresses 
A, while that of c has small effect. This is because the 
reduction in c mainly enhances the dx^-y^ hopping, and 
the suppression of the electron correlation within this or- 
bital has small effect when the intraorbital interaction is 
small. 
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FIG, 12: (Color online) s-wave and d-wave eigenvalues of the 
Eliashberg equation calculated for orbital-dependent interac- 
tions, plotted against (a) n for the model of LaFeAsO, (b) 
ZAs or /lAs for n = 6,1 with the lattice constants fixed at the 
original values for LaFeAsO, (c) a, and (d) c with ^as fixed at 
the original value for LaFeAsO, In (a), the light red (or gray) 
symbol for s-wave at n = 6,2 indicates that the gap is nearly 
nodal, while the red (or solid) symbol at n = 6,3 stands for 
the nodal s-wave. 



Another effect of adopting orbital-dependent interac- 
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LaFeAsO 
orbital dependent interactions 




FIG. 13: (Color online) The gap function for band 4 calcu- 
lated with the orbital-dependent interactions for (a) LaFeAsO 
with n = 6.1 (left) and n — 6.3 (right), and (b) zas ~ 0.658 
(left) and zab ~ 0.6304 (right) with n = 6.1. (c) Xs3333 and 
Xs4444 for 2As = 0.658 and n = 6.1. 



VI. CALCULATION FOR ACTUAL 
MATERIALS 

In this section, we calculate the band structure of ac- 
tual materials other than LaFeAsO, i.e., the phosphate 
and Nd compound, using the experimentally determined 
lattice structure to construct the five-band model. The 
band filling will be fixed mainly at n = 6.1 to make a 
direct comparison with the results for LaFcAsO. The re- 
sults are interpreted in view of the general trend obtained 
in the study of the virtual lattice structures. 



A. LaFePO 

The band structure of the five-band model for LaFePO 
is shown in FiglTH In the case of LaFePO, the lattice 
constants are small compared to LaFeAsO (while closer 
to NdFeAsO below). However, the hopping integrals are 
similar to or larger than those for the virtual structure 
for LaFeAsO with zas = 0.6304. Thus, the main differ- 
ence from LaFeAsO is caused by the height of P. The top 
of the dx^-Y^ band at (tt, tt) is indeed pushed below the 
Fermi levelj^ and this makes the (tt, 7r/2) spin fiuctua- 
tions arising from the /?-/? nesting dominate in Xsna 
shown in FigfTSl 



NdFeAsO 




(0.0,0) (11,0,0) (7l.ll,0)(Jl,7t,7t) (0,0,0) (0,0,JI) 



LaFePO 



tions appears in the symmetry of the gap function. For 
ZAs ~ 0.658 with n — 6.1, we have seen that the gap is 
large at the dx'^-Y^ charactered portions of the Fermi 
surface when we adopt orbital-independent interactions. 
For the orbital-dependent interactions, the absolute value 
of the gap is nearly constant for each of the all pockets 
as shown in the left panel of FigfT^Tb) for band 4. This 
should be again because the magnitude of the gap is re- 
duced at dx^-Y^ portions of the Fermi surface due to the 
reduction in the dx^-Y^ intraorbital interaction. The ef- 
fect of reducing the dx^-Y^ orbital interaction can be 
clearly seen in the comparison between Xs3333 and Xs4444 
depicted in Fig fT5f c). where the two components have 
similar magnitudes for zas = 0.658 in contrast to the re- 
sult for the orbital- independent interactions in Fig[5Ja). 
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FIG. 14: (Color online) The band structure of the five-band 
model of the optimally doped NdFeAsO (upper panels) and 
LaFePO (lower). The Fermi surface at fc^ = for n = 6.1 is 
shown on the right. 



This behavior in the spin fluctuation for LaFePO 
acts to make the d-wave pairing (Fig ll5l right) domi- 
nate for the orbital-independent interactions, while the 
sign-reversing s-wave with nodes intersecting the /? Fermi 
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FIG. 15: (Color online) Upper panels : Xs3333 a-nd X34444 for 
the five-band model for LaFePO with n = 6.1, U = 1.7, U' = 
1.4, and J = J' = 0.15 (orbital-independent interactions). 
Lower panels : s-wave (left) and d-wave (right) gap functions 
for bands 3 and 4 for the same parameter values. 



surface (Fig[Tn]left) is also closely competing. The U de- 
pendence of the eigenvalue of the Eliashberg equation is 
shown in Fig fTCT a). If we adopt the orbital-dependent in- 
teractions introduced in sectior fVCl on the other hand, 
we find that the nodal s-wave slightly dominates over d- 
wave in the entire parameter regime studied, as depicted 
in Fig flCT b) for the eigenvalue of the Eliashberg equation 
against the interaction strength (i.e., the multiplication 
factor / here) . This is expected from the comparison be- 
tween Figin^a) and FigfTWa). where we can see that the 
orbital-dependent interaction tends to favor nodal s-wave 
over d-wave. 

We further find that the s-d competition depends on 
the band filling (not shown), i.e., smaller band fillings 
tend to favor the nodal s-wave. As seen from these re- 
sults the competition between nodal s-wave and d-wave 
pairings in LaFePO is rather subtle, and it is difficult to 
theoretically determine which symmetry actually takes 
place. In either case, however, superconducting gap of 
LaFePO is expected to have nodes intersecting the Fermi 
surface. This is in fact consistent with recent experiments 
on LaFePO that suggest the presence of nodes in the su- 
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FIG. 16: (Color online) s-wave and d-wave eigenvalues of the 
Eliashberg equation for LaFePO plotted against (a) U with 
orbital-independent interactions with U—U' = 2J = 0.3 fixed, 
and (b) the multiplication factor / with orbital-dependent in- 
teractions. The s-wave gap here always has nodes intersecting 
the /3 Fermi surface. 
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B. NdFeAsO 

The band structure of NdFeAsO is shown in FigfT4l 
The low-temperature lattice structure of the optimally 
doped sample (sample 4) in ref 50 is adopted here. Here, 
we have performed an LDA calculation with the plane 
wave basis set using the pseudopotential of Nd obtained 
with the open-core treatment for the / electrons. The 
dx^-Y^ band at (tt, tt) is seen to cross the Fermi level 
even at 71 = 6.1 as expected, and the (tt, 0) spin fluctu- 
ation strongly dominates in Xs4444 shown in the right 
panel of Fig fTTT a). and the fully-gapped, sign-reversing s- 
wave shown in Fig lisr a) strongly dominates over d-wave. 
For the orbital-independent interactions, the eigenvalue 
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of the Eliashberg equation is A = 1.20 for s-wave, which 
is indeed greater than that for LaFeAsO (A — 0.90), but 
not as large as the virtual lattice structure of LaFeAsO 
where the As height is increased to the value of NdFeAsO 
(A — 1.91). The latter property can mainly be attributed 
to the reduction in the lattice constants a and c in Nd- 
FeAsO as compared to those of LaFeAsO. 

(a) optimally doped NdFeAsO 
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Xs4444 



















(b) underdoped NdFeAsO 




FIG. 17: (Color online) Xs3333 (left panels) and X34444 (right) 
in the five-band model for (a) the optimally doped NdFeAsO 
with n = 6.1, and (b) the underdoped NdFeAsO with n = 
6.03. 



When we adopt the orbital-dependent interactions in- 
troduced in sectioi fVCl the enhancement of A from the 
LaFeAsO value is weaker than in the case of orbital- 
independent interactions, as expected from the previous 
discussion. Namely, ALa = 0.70 and ANd = 0.72 for the 
multiplication factor / = 0.42 while ALa = 1-20 and 
ANd — 1-32 for / = 0.45. The magnitude of the gap 
for the s-wave pairing is nearly constant on each of the 
all pockets as shown in Fig fTST b). which is also expected 
from the argument in sectior fV CI 

We have also performed a calculation for NdFeAsO 
with the lattice structure in the underdoped regime (sam- 
ple 1 in reflsoh. Here we take the band filling of n = 6.03 
and adopt orbital- independent interactions. As shown in 
FigflTl the maximum value of XsAAAA (and also Xs, not 
shown) is larger than in the case of the optimally doped 
lattice structure because the nesting is better for the un- 
derdoped case. Nonetheless, the eigenvalue of the Eliash- 
berg equation is found to be smaller, A = 1.03 for the 
s-wave as compared to A — 1.2 for the optimally doped 
sample. This shows that removing the electrons to make 
the 7 Fermi surface more effective does not necessarily fa- 
vor superconductivity. This may be because lowering the 
Fermi level results in the decrease in the dx'^-Y^ density 
of states on the (3 Fermi surface, where the dx^-Y^ band 
forms a Dirac cone. Doping the electrons raises the Fermi 
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FIG. 18: (Color online) The s-wave gap function for band 
3 (left panels) and band 4 (right) for the optimally doped 
NdFeAsO. Orbital-independent (a) and dependent (b) inter- 
actions are adopted. 



level, thereby increasing the dx^-y^ density of states on 
the [3 Fermi surface, and at the same time increasing the 
pnictogen height, which pushes up the dx^-ya band at 
(tt, tt) so as to catch up with the raised Fermi level. The 
increased dx^-y^ density of states upon doping is seen 
in the broad (tt, 0) peak structure in XsAAaa in the opti- 
mally doped case as compared to that in the underdoped 
regime (Fig fTf)) . 

The tendency that electron doping tends to increase 
the pnictogen height should be general because the in- 
crease in the negative charge in the Fe layers suppresses 
the attractive interaction between the positively charged 
iron and the negatively charged pnictogen. Therefore, 
the doping dependence of A shown in Fig[6Ua) with a 
fixed band structure may be too naive in that the effect 
of the 7 Fermi surface monotonically decreases with the 
higher electron doping. As for the doping dependence 
of superconductivity, other than the effect of the change 
in the band structure, there have also been theoretical 
studies that suggest the importance of the "unscreening 
effect" of the Coulomb interaction®^, or the importance 
of the electron correlation^^ 

We have finally examined the effect of pressure on Nd- 
FeAsO. The lattice structure data for NdFeAsOi_j, un- 
der a pressure of 3.8GPa is taken from refisil. Applying 
pressure on NdFeAsO tends to reduce /ias as well as the 
lattice constants. At 3.8GPa, the height is reduced to 
/lAs = I.33A. This results in a suppression of the eigen- 
value of the Eliashberg equation, and the s-wave eigen- 
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value is reduced regardless of the choice of the electron- 
electron interactions, e.g., for the orbital- independent in- 
teractions, A at T = 0.02 reduces (from 1.2) to 0.67. The 
suppression of A is at least qualitatively consistent with 
the experiment i^iS 



VII. DISCUSSIONS 
A. Validity of the five-band model 

One of the most important next steps in the micro- 
scopic study on superconductivity in the iron-based su- 
perconductors should be an examination of the present 
scenario based on the Fermi surface nesting by means of 
self-consistent calculations. Indeed, if we are interested 
in the behavior of the present five-orbital (d-only) model 
with moderate (realistic) size of the Hund coupling^^ or 
the Coulomb interaction a''^'^^ estimated by various ab 
initio methods, we have obviously to take account of 
the self-energy (otherwise we have magnetic ordering at 
rather high temperatures), where the self-energy correc- 
tion generally affects the effect of Fermi surface nesting, 
which is usually overestimated in RPA. 

One possible way to go beyond RPA is employing the 
fluctuation exchange (FLEX) approximation^^. How- 
ever, it has recently been recognized that FLEX does 
not work so successfully for the five-orbital model with 
moderate correlations. Namely, while we naively ex- 
pect that the model should have a strong instability for 
the stripe-type antiferromagnetic ordering for undoped 
LaFeAsO, the spin susceptibility in FLEX has a peak at 
(tTjTt), which corresponds to the checkerboard- type anti- 
ferromagnetic instability. Even in the weakly correlated 
regime, Ikedapi^ had to introduce artificial level shifts for 
and dx'^-y^ to the original LDA band in order to avoid 
a large d^^ Fermi surface. 

These problems seem to come from the fact that the 
self-energy correction generally has a strong orbital de- 
pendence in the five-orbital model rather than from some 
problems in FLEX. Even in the simple Hartree approx- 
imation for the paramagnetic case, the band structure 
and the Fermi surface dramatically change from those in 
LDA due to the Hartree field (~ U{ni)), since the fill- 
ing of each of the five orbitals varies so differently (e.g., 
n - 0.8(0.5) for d^2{d^2_y2)). 

On the other hand, in the dpp model that takes account 
of Fe 3d and As Ap and O 2j»^ for which five Fe 3d orbitals 
are similarly filled, we do not have such problems as far as 
we introduce a gap (the so-called double-counting term 
A) that depends on the difference between the correla- 
tions in Fe "id and As 4p2S,. In fact, this double-counting 
term in the dpp model makes the situation in the five- 
orbital model subtle. For the dpp model, we can safely 
assume that A does not have a serious orbital depen- 
dence. On the other hand, if we translate A in terms 
of the five-orbital model, we have to assume that A has 
a non-trivial orbital dependence, since each Fe 3c? hy- 



bridizes with As Ap differently (e.g., while the hybridiza- 
tion between d^2_y2 and As Ap is strong, those between 
dz2 and As Ap are weak^). This should be one reason 
why Ikcda had to introduce an orbital-dependent level 
shift in his FLEX calculation^'^. 

Therefore, we believe that it is impossible to obtain 
any meaningful results in self-consistent calculation for 
the five-orbital model without considering the orbital- 
dependent double-counting term, while we have still no 
guarantee that the double-counting term can really make 
the five-orbital model mimic the original dpp model. An 
interesting observation is that the situation is in sharp 
contrast with the case of high Tc cuprates. For cuprates, 
aside from the issue of the validity of the single-band 
Hubbard model or the t-J model, we can naively expect 
that the self-consistent solutions of these models can at 
least describe Mott insulator, metallic state with strong 
antiferromagnetic fluctuations, etc. On the other hand, 
for iron pnictides we have to seriously examine whether 
the five-orbital model indeed has a self-consistent solu- 
tion with the stripe-type antiferromagnetic instability. 
This is why we consider this an important future prob- 
lem. 



B. Phase diagram 

In the preceding sections, it has been shown that the 
pnictogen height can act as a switch between high Tc 
nodeless and low Tc nodal superconductivities. We have 
also shown that the increase in the lattice constants is 
unfavorable for superconductivity. These tendencies can 
be incorporated in a schematic phase diagram shown in 
Figim In the upper panel, we take the pnictogen height 
as the horizontal axis, and lattice constants as the ver- 
tical axis. In LnFeAsO, the lattice constants decrease 
monotonically in the chemical trend La— + Nd— + Dy.^^ On 
the other hand, the As height monotonically increases, 
and these effects may cancel with each other to result in 
a nearly constant Tc between Ln=Nd and Dyi^ In the 
lower panel of Fig[T9l we adopt the As-Fe-As bond angle 
a as the horizontal axis to make clear comparison with 
reflsol. in which it was shown that the maximum Tc seems 
to be reached when the pnictogens form a regular tetra- 
hedron. As schematically indicated by a curved arrow, 
the appearance of the maximum Tc may be a consequence 
of the combined effect of the bond angle and the lattice 
constants. Thus, as far as the present theoretical study 
is concerned, the pnictogen height is a better parameter 
than the angle a to draw a phase diagram, since a is 
affected by both the height and the lattice constant a, 
which have opposite effects on Tc- 

To attain higher Tc on the basis of this phase dia- 
gram, it is desirable to have higher position of the pnicto- 
gen while keeping lattice constants not reduced. On the 
other hand, we have to keep in mind that such a varia- 
tion in the lattice parameters also enhances the tendency 
toward magnetism. Experimentally, magnetic ordering 
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FIG. 19: (Color online) Schematic superconducting phase di- 
agram on the hpn-[a,c] (upper) and a-[a,c] (lower) planes. 
The dashed lines are schematic contours of Tc- The curved 
arrow in the lower panel schematically indicates how the lat- 
tice parameters vary as the bond angle a is decreased in refls^ 
The shorthands for materials are the same as in Table IH 



is icnown to occur after the structural phase transition, 
and the tendency toward magnetism and the structural 
phase transition seem to be linked. In this sense, the 
superconductivity will have to compete with the mag- 
netism/structural phase transition in a more severe man- 
ner as the height and/or the lattice constants are in- 
creased. Namely, too much increase in the height and/or 
the lattice constants can be unfavorable for supercon- 
ductivity in that the magnetism/structural phase tran- 
sition, which may take place at high temperatures, can 
dominate over superconductivity. In such a case, apply- 
ing pressure to reduce the lattice constants and/or dop- 
ing carriers can be effective to suppress the magnetism 
or the structural phase transition. This seems to be 
the case for the undoped LaFeAsO'^^ and CaFe2As2r^ 
where structural phase transition and magnetism take 
place at ambient pressure, while applying pressure seems 
to remove the phase transition to result in superconduc- 
tivity. A better understanding of the competition be- 
tween superconductivity and magnetism should require 
further understanding of the magnetic state sitting next 
to the superconducting state. Namely, the mutual rela- 
tion among the so-called "spin density wave" state ob- 



served experimentally^!^, the magnetic state obtained 
in first principles calculations^i2£i2ii2^, and the spin fluc- 
tuation obtained in the downfolded five-band models has 
to be made clearer. 

Another point that should be kept in mind is the ma- 
terial dependence of the electron-electron interactions, 
which is not considered in the present phase diagram. In 
particular, the difference in the lattice structure in the 
11 systems such as FeSe^ and in the hole doped 122 sys- 
tems such as BaFe2As2^ can affect the effective electron- 
electron interaction within the FeAs planes. Also, the 
screening effect of the / orbitals may be different between 
LaFeAsO and NdFeAsO. On the / electrons, the present 
analysis adopts open-core treatment for the / electrons 
of Nd, whose validity may have to be examined. We 
cannot completely rule out the possibility that the hy- 
bridization between Nd/ and Fed electrons can affect the 
electronic states and thus the superconductivity. In these 
senses, the actual phase diagram for the entire family of 
the iron-based superconductors should have more axes 
than presented here. 



VIII. CONCLUSION 

In the present study, we have investigated how the lat- 
tice structure affects the spin-fluctuation mediated super- 
conductivity in iron pnictides. The obtained picture is 
that the gap function and Tc are determined by the com- 
petition or cooperation of the multiple spin fluctuation 
modes arising from a-/3. Pi- 02, and 7-/3 Fermi surface 
nestings, which depends on the materials, band filling, 
and/or pressure. In particular, the competition between 
and 7-/3 nestings within portions of the Fermi sur- 
face having strong d^^-y^ character governs the form 
of the superconducting gap as well as the strength of 
the superconducting instability. The relevance of the 
7 (quasi) Fermi surface is determined by the pnicto- 
gen height, and consequently, the pnictogen height plays 
the role of a switch between high-Tc, nodeless and low- 
Tc, nodal pairings, which may give the answer to the 
question of why the form of the superconducting gap 
as well as Tc are vastly different between LaFeAsO and 
LaFePOi ^^'^^ An intriguing observation is, since d-wave 
and the nodal s-wave tend to be closely competing for 
low pnictogen heights, there is a possibility of exotic 
pairing such as s -I- ic?.^^ The lattice constant also af- 
fects superconductivity in the manner that the reduction 
in a (c) mainly suppresses the electron correlation within 
dxz/dyz {dx^-Y"^) orbitals and thus degrades supercon- 
ductivity. A schematic phase diagram has been obtained 
by combining the effect of the pnictogen height and the 
lattice constants. 

The low Tc in cases where the /3i-/32 nesting dominates 
over 7-/3 nesting can be naturally understood in the sense 
that there is a kind of frustration between P1-P2, a-/3i, 
and a-/32 nestings in determining the form of the gap be- 
cause the sign of the gap has to be changed across each 
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of the multiple nesting vectors (see Fig|4]). Conversely, 
the high in the case where the 7-/3 nesting, along with 
a-(3, dominates over /3i-/32 is natural in that the unfrus- 
trated gap fully opens on all five disconnected pieces of 
the Fermi surface, ai, a2, Pi, P2, and 7. In this sense, 
high Tc in iron pnictides can be understood as a realiza- 
tion of the theoretical proposal that one can look for high 
Tc superconductors in systems with disconnected Fermi 
surfaces 

The form of the s-wave gap is shown to be nonuniver- 
sal, even when the gap is fully open, and its variation 
along the Fermi surface strongly depends on the band 
structure (i.e. the lattice structure), the band filling, and 
the electron-electron interactions. When the s-wave gap 
varies significantly along the Fermi surface, it is also ex- 
pected to be affected by the presence of the impurities, as 
pointed out by Mishra et al.^^ In this sense, the form of 
the gap should experimentally be determined by a com- 
bination of multiple experiments on the same material, 
desirably on the same sample. From this viewpoint, the 
discrepancy between the NMR experiments for LaFeAsO 
and other experiments suggesting nearly isotropic gap 
may be a consequence of the nonuniversality of the su- 
perconducting gap, especially because LaFeAsO hes close 
to the nodeless/nodal boundary. In fact, anisotropic 
sibwave pairing where the gap varies strongly on the /? 
Fermi surface has been proposed to explain the T"^ decay 
in the NMR experiment This view may also give some 
clue as to why some of the tunneling spectroscopy mea- 
surements exhibit zero bias conductivity peak,^^ which is 
an indication of unconventional sign reversing pairing^iS^ 
while others do noti^Sl It is worth noting that recent the- 
oretical studies show that it is unlikely to observe the zero 
bias conductivity peak for the fully gapped sign reversing 
s-wave pairingJ£2ii2^ 



In the present study, we have focused on the LnFeAsO 
(1111) systems. In the 11 systems such as FeSe, ex- 
periments under pressure also suggest strong structure 
dependence of superconductivity ;^°^i^°^'^°^i^°^ However, 
there may be some discrepancies from the 1111 systems 
regarding the lattice structure dependence since the LnO 
layer is not present. Also, the extremely high position of 
the chalcogen atom in 11 systems (1.47A in FeSe^"^ and 
1.76 A in FeTel°?) may affect the lattice structure depen- 
dence. Our study focusing on the 11 system is now un- 
derway. Similarly, the hole-doped 122 systems are also 
expected to have some discrepancies with the electron- 
doped 1111 systems, which also deserves future study. 
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